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Abstract. The Hipparcos Intermediate Astrometric Data for carbon stars iiave been reprocessed using an algorithm 
which provides an objective criterion for rejecting anomalous data points and constrains the parallax to be positive. 
New parallax solutions have been derived for 317 cool carbon stars, mostly of types R and N. In this paper we 
discuss the results for the R stars. The most important result is that the early R stars (i.e., RO - R3) have absolute 
magnitudes and V — K colors locating them among red clump giants in the Hertzsprung-Russell diagram. The 
average absolute magnitude Mr for early R-type stars (with V — K < 4) has been derived from a Monte-Carlo 
simulation implicitly incorporating all possible biases. It appears that the simulated magnitude distribution for 
a population with a true Gaussian distribution of mean Mk = —2.0 and intrinsic standard deviation 1.0 mag 
provides a satisfactory match to the observed distribution. These val ues are cons istent with the average absolute 



magnitude Mk = —1.6 for clump red giants in the solar neighborhood (Alves 2000). Further, early R-type stars are 
non-variable, and their infrared photometric properties show that they are not undergoing mass loss, properties 
similar to those of the red clump giants. 

Stars with subtypes R4 - R9 tend to be cooler and have similar luminosity to the N-type carbon stars, as confirmed 
by their position in the {J — H,H — K) color-color diagram. 

The sample of early R-type stars selected from the Hipparcos Catalogue appears to be approximately complete 
to magnitude Kq ^ 7, translating into a completeness distance of 600 pc if all R stars had Mk ~ —2 (400 pc if 
Mk = — 1). With about 30 early R-type stars in that volume, they comprise about 0.04% (0.14% for Mk = — 1) of 
the red clump stars in the solar neighborhood. Identification with the red clump locates these stars at the helium 
core burning stage of stellar evolution, while the N stars are on the asymptotic giant branch, where helium shell 
burning occurs. The present analysis suggests that for a small fraction of the helium core burning stars (far lower 
than the fraction of helium shell- burning stars), carbon produced in the interior is mixed to the atmosphere in 
sufficient quantities to form a carbon star. 
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1. Introduction 

The Harvard spectroscopic classification of stars for the 
Henry Draper Catalogue recognized the first examples of 
chemically-peculiar stars, the carbon and S giants. Carbon 
stars are immediately recognizable by the presence of ab- 
sorption bands of C2 at 438 3 A and 4737 A and the 4216 
A band of CN ( Shane 192^), and are chem ically charac- 
terized by having C/0 > 1. dickering (1896, 1908) showed 
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that the carbon stars further subdivide into two spectral 
classes, the 'warm' R stars and the 'cool' N stars, and in 
the intervening years it has become increasingly apparent 
that these stars do not represent a continuum of prop- 
erties but rather belong to two distinct classes. Further, 
Keenan (1942| ) noted that some members of the R spec- 
tral class were metal-poor and had high radial velocities 
and proper motions. These are now separately classified 
as CH stars. Several low-concentration globular clusters 
contain CH stars ( [Bond 197^ ; [McClure 1984bD ; they are 
generally considered to be halo carbon stars. The most 
recent list of CH stars is by Bartkevicius (1996| ). Several 
additional types of carbon star have been recognized, in- 
cluding R CrB stars and dwarf carbon stars. 
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The revised classification system of Keenan (1993) is 
an attempt to provide consistent temperature, luminosity 
and abundance sequences among the different types of car- 
bon stars. Three main classes are defined in this scheme: 
C-Rti, C-Nn, and C-Hrt, corresponding to the old R, N 
and CH types respectively, with n being a homogeneous 
temperature index for the three sequences, corresponding 
to the G-K-M temperature sequence in "oxygen" stars. 
Unfortunately, spectral types based on this scheme are 
available for only a small fraction of carbon stars. 

Since C/0 < 1 everywhere in the interstellar medium, 
in the Sun, and in most stars, the excess carbon in carbon 
stars must be the result of stellar nucleosynthesis in the 
star itself or in a binary companion. The many varieties of 
carbon stars are thus stars in the late stages of evolution 
or wit h evolved binary companions, and are preferentially 



stars have the same absolute magnitude and color, and are 
located in the red clump region of the HR diagram. The 
Hipparcos parallaxes also allow us to derive the spatial 
density of the red clump stars and the fraction of these 
which are R stars; we show that they are very rare, far 
less than 1% of the red clump stars (Sect. 4). The carbon 
enhancements in R stars may therefore originate in the 
helium-burning interiors, perhaps as the result of a mixing 
in a rapidly rotating star spun up by the accretion of its 
companion (see Sect. ^). Sect. 2 discusses the astrometric, 
photometric, spectroscopic and variability data on which 
the results presented in this paper are based. The absolute 
magnitudes and colors are examined in Sect. 3, the space 
densities in Sect. 4, and the implications and conclusions 
are discussed in Sect. 5. The details of the re-reduction of 



the Hipparcos data will be discussed by Pourbaix et al 



formed in systems whose intrinsic metallicities (including (2001a, Paper III) and the results for N-type carbon stars 



the oxygen abundance) are low. 

The N stars are understood in a general sense. They 
are very luminous, cool stars with shell hydrogen and he- 
lium burning and with Sa-produced ^^C mixed to the sur- 
face by deep convection ( Iben fc Renzini 1983| ), although 
the details of how that mixing occurs are far from well un- 
derstood. Other stellar types with enhanced carbon abun- 
dance, like CH and barium stars, are observed to be mem- 



by pourbaix et al. (2001b| ,Paper IV) 



bers of binary systems (see the discussions by McClure 
1984a[] [i997a ; Jorissen et al. 1998), and their peculiar 



abundances can be attributed to mass transfer. The R 
CrB stars are variable hydrogen-deficient stars, and have 
shed most of their envelope. The dwarf carbon stars ( Dahr| 
et al. 1 377 ; Green & Margon 1994) and the R stars remain 



a mystery. Like the N stars, the R stars are red giants, but 
their average luminosity, determ i ned by stat i stical paral - 
lax methods ( Vandervoort 1958 ; Scalo 1976 ; Dean 1976 ) 
is too low for shell helium burning. Furthe r, their oxyge n 
abundances are similar to that of the Sun ( Dominy 1984 ), 
they do not exhibit s-process overabundances, and no R 



star has been found so far in a binary system (McClure 
1997b^ a statistically unlikely and surprising result. It 
has therefore been suggested that R stars might actually 
be coalesced binary systems. The remaining hypothesis is 
that these stars are helium core-burning stars or post he- 
lium core burning stars, with carbon produced by helium 
burning somehow introduced into their atmospheres, per- 
haps at the helium core flash (Scalo 1976; Dominy 1984). 
There is little direct evidence for this hypothesis, though, 
because of lack of information about the intrinsic luminos- 
ity function of the R stars or of their spatial density. This 
paper sets out to derive the luminosities of the R stars 
to attempt to discover their location on the Hertzsprung- 
Russell diagram and establish their evolutionary status. 



The Hipparcos project (ESA 1997) made it possible 
for the first time to make direct estimates of the dis- 



tances to giant carbon stars. However, as noted by [Knapp 
(200l|)] the parallax errors for these stars are often far 
larger than expected. In this paper, we use new parallaxes 
derived from a re-reduction of the Hipparcos Intermediate 
Astrometric Data to show that essentially all the (early) R 



2. The data 

2.1. Overview 

The basic data for R stars, including the new paral- 
lax data, are collected in Tables |l] & |. Stars with en- 
hanced carbon abundances come in many spectral types: 
the N stars, the R stars, carbon white dwarfs, carbon- 
rich Wolf-Rayet stars, carbon cepheids, R CrB stars, bar- 
ium stars, CH stars, and dwarf carbon stars, among oth- 
ers. As pointed out above, their chemical peculiarities are 
due to nucleosynthesis within the star itself, or its com- 
panion, and the task is to identify chemical peculiarities 
with particular evolutionary paths. This is carried out in 
the first instance by associating like objects, classifying 
them by spectral type, luminosity, color, and variability, 
and is made difficult both by the wide variety of stars 
with carbon enhancements, and by the often uncertain 
observational information. A given star can be classified 
as several different spectral types in several different pa- 
pers, causing confusion in attempts to evaluate the com- 
mon properties of given classes of stars (see the discus- 
sions by [McClure 19974 |b|; [Wallerstein fc Knapp 1998[ ). 



Examples in clude HIP 12028 (H P 161115) which is classi- 
fied as R by [Vandervoort (1958[), as CH by [Keenan (1993|) 



and as J by [Ohnaka fc Tsuji (1999), while HIP 8575 



(BD-H02° 3336) is classified as R2 by [Vandervoort (1958[) 
as CH by Yamashita (1975| ), and as N by Barnbaum et al. 
[(1996[) . Thus it is far from straightforward to decide, from 
the available data, whether the presence of a star in a given 
variability class with a spectral type atypical of that class 
demonstrates a range of spectral types in that variability 
class or is due to inadequate data. The infrared color-color 
diagram {J — H, H — K) (Fig. 0) may nevertheless be of 
some help in the process of identifying misclassified stars, 
as the different types of carbon stars occupy specific loca- 
tions in that diagram (see Sect. 2.2). 

Of the many varieties of carbon star, the cool lumi- 
nous carbon stars are of interest because they are ex- 
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Table 1. Early R stars: Astrometric, photometric and spectroscopic data. The rederived parallax Wo and its lower and 
upper bounds (respectively vjl and tjJjj] see text) are given in milliarcseconds (mas). The R stars have been separated 
into early and late types, according to the V — K color index being smaller or larger than 4 (see Fig. ||). 



HIP 


GCGCS 








V 


K 


J-H 


H -K 


Ref 


Sp. 


Ref Rem. 






(mas) 


(mas) 


(mas) 










JHK 


type 


SpT 


2700 


80 


1.32 


2.37 


4.24 


9.62 


6.94 


0.51 


0.16 


D86 


R2 


B96 


5809 


196 


0.24 


1.09 


5.05 


10.07 


7.22 


0.77 


0.03 


C79 


RO 


HD 


7816 


256 


— 


_ 


_ 


11.11 


7.90 


1.02 


0.17 


D86 


RO 


S69 


11508 


357 


0.25 


1.39 


3.78 


9.57 


5.92 


_ 


_ 


D86 


RO 


V58 


12028 


378 


3.52 


4.57 


5.93 


8.15 


5.38 


0.62 


0.15 


F92 


R2 


V58 


17933 


576 


0.54 


1.11 


2.29 


8.29 


4.65 


0.73 


0.19 


C81 


R3 


HD 


18564 


588 


0.37 


1.49 


4.51 


10.26 


6.66 


0.41 


0.23 


D86 


R3 


S28 


18696 


594 


_ 


_ 


_ 


10.45 


7.23 


0.58 


0.15 


2M 


R2 


S28 


19050 


610 


_ 


_ 


_ 


10.71 


8.10 


0.48 


0.16 


2M 


RO 


V58 


19269 


639 


0.72 


1.74 


4.20 


10.62 


6.97 


0.56 


0.46 


C79 


R3 


HD 


26927 


1035 


0.17 


1.08 


6.57 


10.02 


7.62 


0.5 


0.09 


D86 


RO 


V58 


28172 


1110 


0.27 


1.33 


4.40 


10.35 


7.62 


0.53 


-0.11 


D86 


R2 


S28 


31725 


1380 


0.46 


1.10 


2.66 


9.41 


6.46 


0.54 


0.21 


C79 


RO 


HD 


33042 


1460 








10.56 


_ 


_ 


_ 


_ 


R 


S72 


35681 


1622 


0.37 


1.05 


2.99 


8.52 


6.16 


0.4 


0.18 


2M 


RO 


V58 * 


36086 


1703 


0.13 


1.08 


4.48 


9.22 


5.81 


0.61 


0.34 


D86 


R3 


S28 


39118 


1981 


0.06 


1.32 


7.23 


10.41 


7.23 


0.79 


0.25 


2M 


R2 


V58 


40374 


2126 




_ 


_ 


11.70 


8.38 


0.6 


0.22 


2M 


R 


GCGCS 


44172 


2396 


0.81 


2.08 


5.34 


9.44 


6.84 


_ 


_ 


D86 


R2 


B96 


48329 


2626 


0.35 


1.22 


4.23 


9.90 


6.25 


0.7 


0.22 


D86 


R3 


S28 


50994 


2759 


0.35 


1.05 


3.14 


9.56 


_ 


_ 


_ 


_ 


RO 


HD 


53354 


2892 


_ 


_ 


_ 


10.54 


7.62 


0.61 


0.17 


2M 


R2 


V58 


53522 


2900 


0.63 


1.81 


5.15 


10.14 


_ 


_ 


_ 


_ 


RO 


V58 


53810 


2925 


_ 


_ 




8.35 


_ 


_ 


_ 


_ 


K5R 


S44 


53832 


2919 


0.96 


2.08 


4.52 


10.11 


_ 


_ 


_ 


_ 


RO 


V58 


56405 


3058 


0.12 


1.19 


6.52 


10.27 


_ 


_ 


_ 


_ 


R2 


S28 


58786 


3156 


0.26 


1.06 


4.20 


10.27 


7.51 


0.52 


0.15 


2M 


R2 


V58 


62944 




7.41 


8.22 


9.12 


6.91 


4.16 


0.55 


0.1 


E78 


R3 


B96 


63955 


3335 


1.44 


2.26 


3.54 


8.50 


5.97 


0.46 


0.18 


2M 


R3 


B96 


65320 


3379 


1.84 


2.94 


4.71 


9.66 


_ 


_ 


_ 


_ 


RO 


S28 


66317 


3409 


— 




— 


12.39 


_ 


_ 


_ 


_ 


R 


S72 


68543 


3469 


0.21 


1.02 


5.07 


9.48 


6.67 


0.52 


0.19 


D86 


R2 


S28 


69089 




0.53 


1.17 


2.59 


8.68 


6.40 


0.46 


0.12 


2M 


R2 


B96 


71464 


3513 


_ 


_ 


_ 


12.17 


_ 


_ 


_ 


_ 


R 


S72 


73955 


3591 


0.42 


1.17 


3.28 


9.99 


_ 


_ 


_ 


_ 


RO 


S28 


74179 


3562 


0.28 


1.16 


4.83 


9.72 


6.50 


0.5 


0.5 


F97 


R3 


HD * 


74826 


3586 


0.52 


1.55 


4.64 


9.78 










RO 


S28 


75691 


3614 








12.88 


11.12 


0.4 


0.06 


2M 


R5 


V58 * 


75745 


3605 








11.79 










R 


S72 


80365 


3687 








11.18 


8.16 


0.23 


0.5 


F73 


R 


M51 * 


80769 


3704 


0.37 


1.37 


5.11 


10.40 


7.50 






M65 


R2 


V58 


82184 


3735 


0.76 


1.74 


3.99 


9.10 


6.46 


0.57 


0.12 


2M 


R2 


B96 


84266 


3795 


2.30 


2.89 


3.62 


7.60 


5.11 


0.3 


0.23 


N81 


R2 


B96 


85117 


3816 


1.88 


3.33 


5.91 


9.57 










R3 


C21 


86927 


3879 


3.77 


4.78 


6.05 


8.71 


6.14 


0.53 


0.14 


D86 


RO 


S28 


87603 


3912 


0.30 


1.05 


3.66 


10.72 


7.92 


0.55 


0.18 


2M 


R4 


S44 


88584 


3938 


2.07 


3.49 


5.88 


9.82 










R3 


S44 * 


89239 


3973 


4.56 


6.37 


8.91 


10.82 


7.41 


0.61 


0.24 


2M 


R5 


V58 



pelling carbon-enhanced mass into the interstellar me- 
dium. Knapp (2001 ) identified about 320 carbon stars in 
the Hipparcos catalogue. An attempt was made to be com- 
plete (as far as the members of the Hipparcos catalogue 
are concerned) for the luminous cool carbon stars, i.e., the 
R and N stars (because these are the stars which are, or 



may be, losing mass). Although the Hipparcos catalogue 
itself is based only partly on a complete input catalogue 
(the Input Catalogue contained both complete samples of 
stars - designed as 'Hipparcos survey starsj^ - and addi- 



1 Only 3 such stars (HIP 12028, HIP 62944 and HIP 84266 
discussed herein) are present in the sample of R stars 
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HIP 


GCGCS 


■COL 


■CUo 


TaJu 


V 


K 


J-H 


H ~K 


Ref 


Sp. 


Ref 


Rem. 






(mas) 


(mas) 


(mas) 










JHK 


type 


SpT 




89739 


3982 


2.53 


4.61 


8.43 


10.67 


8.00 


0.5 


0.5 


F97 


RO 


C21 




90199 


4002 


- 


- 


- 


11.91 


- 


- 


- 


- 


RO 


S44 




90694 


4021 


0.15 


1.19 


5.56 


9.88 


6.30 


0.67 


0.23 


C79 


R5 


S44 




92207 


4098 


0.26 


1.62 


6.91 


10.00 


6.5 


1 


1 


F97 


RO 


HD 




93987 


4181 


0.33 


2.09 


13.40 


11.08 


5.89 


0.54 


0.52 


2M 


R2 


S44 




95392 


4485 


0.39 


1.06 


2.86 


9.66 










R2 


S28 




98057 


4560 


0.24 


1.12 


5.23 


9.62 










RO 


S28 




99725 


4784 


0.48 


1.10 


2.49 


9.75 


6.05 


0.7 


0.15 


U97 


RO 


V58 




102726 


4972 








10.64 










Rl 


S44 




105212 


5227 


0.98 


2.31 


5.44 


9.71 










R2 


V58 




105241 


5230 


1.56 


3.11 


6.19 


9.81 


6.97 


0.5 


0.15 


D86 


R2 


V58 




107349 


5408 


0.20 


1.04 


5.42 


10.16 


5.51 


0.91 


0.36 


C79 


RO 


HD 




114452 


5822 


0.72 


1.70 


4.00 


9.99 


7.03 


0.47 


0.2 


D86 


R2 


S28 




114509 


5823 


1.92 


3.24 


5.47 


9.30 


6.69 


0.96 


0.11 


D86 


RO 


S44 




117467 


5937 


0.98 


1.97 


3.95 


8.52 


5.42 


0.61 


0.21 


C79 


R3 


B96 




Notes 



HIP 35681: RU Cam, type CWa, period 22 d; HIP 74179: S Aps, type R CrB; HIP 75691: NSV 7110; HIP 80365: RT Nor, type 
R CrB; HIP 88584: W CrA, type SRb, period 125 d; HIP 89739: RS Tel, type R CrB; HIP 92207: V CrA, type R CrB; HIP 
93987: SV Sge, type R CrB 

References to spectral types and magnitudes: 

HD: H enry Draper Catalogue; GC GCS: Stephenson (1989|); 2M: Two Micron All-Sky Surve y (|Skrutskic 19971); B96: [Barnbaum 
et al. 
Elias 



(|i996|) 



F73: 



C21: [Cannon (19211); C79 : |Catchpolc ct al. (1979|); C81 



Mendoza V & Johnson (196! 



Cohen et al. (198l|); D86 : [Dominy 



ct al. 



(19861 ); E78: 



Feast & Gl ass (1973|); F92: [Feast fc Whitclock (1992D; F 97: [Feast et al. (19971) ; M5 1: [Mayall (195l[); M65 
N81: pSfoguchi et al. (1981^ ; S28: [Shane (1928|); S44 : panford (1944D ; S69: plettebak et al. (1969) ) 
S72: [Stock & Wroblewski (197^); U97: [UUa et al. (19971); V58: [Vandervoort (195g|) 



tional samples of stars of particular interest), it turns out 
to be reasonably complete for carbon stars to Hipparcos 
magnitude Hp - 9.5 ([Knapp 200l|) or to A' - 7 (Fig. |). 



region as N-type carbon stars, which is not surprising as 
most - if not all - late R-type carbon stars may actually 
be N stars (see the discussion in Sects. 2.3 and 2.4). 



2.2. Photometry 



2.3. Variability 



We use the V magnitude given in the Hipparcos catalogue 
(field H5). We acquired J, H and K photometry from the 
literature, as listed in Tables |l| & § along with the bib- 
liographic reference. This was found with the aid of the 
electronic version of the literature compilation catalogue 
of Gezari et al. (199^ ). In addition, we extracted Ks pho- 



tometry from the Two Micron All Sky Survey (2MASS) 
on-line data release ( ^krutskie 1997 ). In many cases, the 
stars are bright enough that they are saturated in the 
2MASS data and we use older photometry. There are no 
systematic offsets evident among the various photometric 
measures for a given star, and the measurements (includ- 
ing the narrower 2MASS Kg band magnitudes) agree to 
better than about 0.2 mag in almost all cases. We were 
able to find K magnitudes for all but 24 stars. Although 
these observations were made at diff erent epochs, R stars 
are not in general variable (Sect. 2.3) and so there is little 
or no uncertainty introduced into the V or K magnitudes 
by variability. 

Fig. presents the {J — H,H — K) color-color diagram. 
It is immediately apparent that R stars occupy a very 
specific region of that diagram, which is also populated 
by CH stars. A few R stars appear to occupy the same 



N stars are in general long-period variables: irregular (Lb), 
semi-regular (SR or SRb) o r Mir a (Mira or SRa) - see 
Kerschbaum & Hron (1994) and Wallerstcin fc Knapp| 



1998 ) for a discussion of the SRa and SRb stars. R stars 
generally do not appear in variable star catalogues. Of 
the 83 R stars in our sample, 18 appear in the General 
Catalogue of Variable Stars (Kholopov et al. 1998), but 
most of these are late R stars, as listed in Table y. Most 
of the R stars of variability types Mira and SR are also 
losing mass, as indicated by their excess 12 /zm emission 
(Table ^), and are therefore probably closely related to N 
stars. 



2.4. Spectral Types 



There are two major studies of R stars. Vandervoort 



(1958) observed a sample of 98 R stars selected from an 



objective prism survey. McClure (1997b) discusses a sam- 
ple of 22 R st ar s. We have also used sp ectral types from 
Keenan (1993| ), Barnbaum et al. (1996 ), and other refer- 
ences. In toto, we have a sample of 83 R stars. The data 
are listed in Tables fc ||. Where available, we have listed 
the spectral type from RO to R8. 
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Table 2. Late R stars: Astrometric, photometric and spectroscopic data. Columns have the same meaning as in 
Table |l|. Additional columns provide the variability name, type and period (from GCVS), and the K — [12] (= K + 
2.51ogF12/28.3) color index (where F12 is the flux density in the IRAS 12 ^m band expressed in Jy). The K - [12] 
index is a mass-loss indicator, non-mass-losing stars having K — [12] 0.7 (e.g. Jorissen & Knapp 1998). 



HIP 


GCGCS 


(mas) 


tJ7o 

(mas^ 


(mas) 


V 


K 


J-H 


H - 


K Ref Sp. 
JHK type 


Ref 
SpT 


GCVS 


Var 


P 
(d) 


K - 


[12] Rem 


35810 


1686 


0.38 


1.26 


4.18 


9.01 










R8 


S28 


V758 Mon 










36623 


1737 


0.38 


1.36 


3.58 


8.02 


2.95 


0.72 


0.38 


N81 


R9 


S28 
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Notes 



1995) 



HIP 36623: also classified as symbiotic variable; HIP 108205: Li star (Boffin et al 
References to spectral types an d magnitudes: 

HD: Henry Draper Cat alogue: IRC: |Ncugcbaucr fc Leighton (1969|): 2M: Two M icron All-Sky Survey { J^krutskio 19971): B96 



Barnbaum et al. (19961): C79: ICatchpolo et al. (1979|): D86: iDominy et al. (198^): K93: j^eenan (199^ ); K 94: jferschbaum & 
Hron (1994); L 44: [Lee et al. (1944|); M65: Mendoza V fc Johnson (19651): N81: [Noguchi et al. (198l|): S28 : ^hane (1928] ); S44: 



Sanford (1944| ); S72: [Stock fc Wroblewski (1972| ); V58: [Vandervoort (1958| ); WOOa: [Whitelock et al. (2000D 



Many of the remaining 317 carbon stars are spectro- 
scopically classified only as carbon stars, and some of these 
may be R stars also. However, as the {J ~ H, H — K) dia- 
gram shows (Fig. not many non-R stars, apart from the 
CH stars, occupy the location of R stars in that diagram. 
CH stars must be treated separately from the R stars, as 
all CH stars are believed to be binary systems (McClure 
1997a||| and must therefore be processed with a binary 
astrometric model, as presented in Paper I (Pourbaix & 
Jorisse fi 2000| ). 



Further, since there is often disagreement in the vari- 
ous references as to the spectral type of a given star (ex- 
amples were given in Sect. 2T), the 'spectral purity' of the 
sample of R stars collected in Tables ^ & || is best assessed 
by looking at the {J — H,H ~ K) diagram. 

The main source of confusion appears to be the fact 
that several late-type R stars from Table || are located in 
the region occupied by N stars. This is not totally sur- 
prising, as late R stars bear similarities with, or might 
even be identical to, N stars. In the following analysis, it 
will therefore be necessary to treat separately R stars with 
F-if < 0.3 [or {V-K)o < 4; see Fig. || from those with 
H — K > 0.3. Finally, the few R stars also classified as 
R CrB variables should be treated separately as well. 



To stress the fact that the early- and late-type R stars 
in the present sample must be considered separately, they 
have been split into two separate lists (Tables |l| & || re- 
spectively). The assignment of a star as an early- or late- 
type R star is mainly based on the {V — K)q index, with 
the threshold set at 4. The {V — K)q index has preference 
over the R spectral subtype in case of conflicting classifica- 
tions. When the {V — K)q index is not available, a criterion 
based on the spectral type is used (with the threshold set 
at R3.5). It will be shown below that, with this choice for 
the thresholds, the two subsamples of R stars present ho- 
mogeneous properties with respect to the R spectral sub- 
type, to mass loss and to variability type, most late- type R 
stars being long-period variable stars, with circumstellar 
dust emission indicative of mass loss. The K — [12] color 
index (with K - [12] = K + 2.5 logF12/28.3, F12 being 
the flux density in the IRAS 12 ^m band, expressed in 
Jy) is considered as a good indicator of (dust) mass loss, 
according to e.g.. Fig. 21 of [Jorissen fc Knapp (1998 ). 
Non-mass-losing stars have K — [12] ^ 0.7. 

Late- type R stars are thus likely to belong to the mass- 
losing AGB, and are therefore not considered further in 
the following discussion, which focuses on R stars with 
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Fig. 1. Top panel: The [J ~ H,H — K) color-color dia- 
gram for the carbon stars from the Hipparcos catalogue 
with J, H and K collected from the literature (as listed 
in Tables |^ & || for R stars). Symbols are as listed in 
the figure, with 'Hd' and 'eHe' standing for 'hydrogen- 
deficient' and 'extreme-helium' carbon stars. Hydrogen- 
deficient carbon stars are similar to R CrB carbon stars, 
except that they do not undergo obscuration episodes 
( Feast et al. 1997| ). Extreme-helium carbon stars are sim- 
ilar to, but hotter than, R CrB carbon stars (Cottrell & 
Lawsoi |i 199^ ). 

Bottom panel: The {{V — K)o,H — K) diagram for the 
sample of R stars considered in the present study 



2.5. Parallaxes 



The Hipparcos catalogue (|ESA 19971) gives parallax values 
at the 2ct level or better [i.e. vo/eiw) > 2, where e{ix7) is 
the quoted parallax uncertainty] for 17 R stars, of which 
15 are of types RO - R2. The resulting absolute visual 
magnitudes for the early R stars cover a range of My = 
—2.8 to H-4.1 (Knapp 2001). However, the parallax errors 
for many of the stars are larger, sometimes much larger, 
than the errors expected for the magnitudes and ecliptic 
coordinates of the stars. Of the nearly 320 carbon stars in 
the sample of [Knapp (2001] ), only 25% (and only 18% of 
the R stars) have "2(t" parallaxes. 

The Hipparcos Intermediate Astrometric Data (IAD; 
van Leeuwen fc Evans 199^ ) for the entire sample of car- 
bon stars were therefore reprocessed as described and 
Papers I & HI. Therefore, only a brief outline of the 
method is given here. 

The position, distance modulus and proper motion are 
the five parameters in the model fit to be adjusted by 
minimization of the IAD abscissa residuals. Included in 
this fitting procedure is the rejection of outliers, defined 



as having a fit residual greater than three times the aver- 
age of the residuals. In practice, two or three observations 
at most were rejected. Fitting the distance modulus en- 
sures that the derived parallax m is always positive, as 
described in Paper I, but at the expense of asymmetric 
error bars and of a bias on w. This bias is of any impor- 
tance only when e(w)/'Dj ^ 1, as shown in Paper HI. The 
distance modulus is better behaved, and must be used for 
computing average absolute magnitudes, as further dis- 
cussed in Sect. ^. 

Another major difference with respect to the reduction 
scheme applied by the Hipparcos consortia resides in the 
bootstrap algorithm used to reject those data points which 
yield inconsistent FAST and NDAC solutions, as described 
in detail in Paper III. 

Monte-Carlo simulations have shown that this repro- 
cessing scheme of the IAD fails when the noise on the 
IAD is too large with respect to the true parallax, in that 
it yields absolute magnitudes that are too faint by as much 
as 5 magnitudes (i.e., parallaxes come out too large). This 
is not really unexpected, as the bootstrap method will 
remove 'discrepant' observations until a solution passing 
through the remaining IAD is finally found. The parallax 
of that solution will clearly be of the order of the noise 
on the abscissa residuals, though the spurious nature of 
that solution may be noticed by the large uncertainty on 
the parallax. For the sample of R stars considered in the 
present paper, such a bias on the absolute magnitude be- 
comes clearly apparent for stars with K > 8, although the 
parallax of somewhat brighter stars is a lready biased as 
well, as seen in Fig. |2| (see also Sect. 3T). In an attempt 
to keep as many unbiased stars as possible, the follow- 
ing criterion has been adopted to reject biased parallaxes: 
K > 8 (or > 10.5 if no if magnitude is available) or the 
star has a simulated absolute magnitude that lies further 
than 4(7 away from the true absolute magnitude of —2 on 
Fig. 1^, a being the standard deviation of the simulated 
absolute magnitudes for stars with K < 7 (i.e., which are 
not biased). Therefore, Tables |l| & || hst only those par- 
allaxes that may be considered as unbiased according to 
the above criteria. 

The results are listed in Tables |l| & ^, which pro- 
vides the HIP number from the Hipparcos catalogue in 
the first column. The star number in th e General Catalog 
of Cool Galactic Carbon Stars (GCGCS; Stephenson 1901 
is listed in column 2. Columns 3 to 5 provide the paral- 
lax wq derived from the minimization on the distance 
modulus, the corresponding lower and upper bounds {vjl 
and wu) computed from Eq. (9) of Paper I. The follow- 
ing columns list the apparent V magnitude (field H5 of 
the Hipparcos catalogue), the K magnitude, the J ~ H 
and H — K color indices, their reference, and finally the 
spectral type and its reference. 
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Fig. 2. The absolute magnitudes M k re trieved by pro- 
cessing the IAD as indicated in Sect. for the sample 
stars assuming that they all have a true absolute magni- 
tude M*^ — —2.0 (see also Sect. 3^). It is clearly seen 
that a strong bias on Mk appears for stars with K ^ 8 
that have noisy IAD. Stars with negative parallaxes in the 
Hipparcos catalogue are identified by filled squares 

3. The absolute magnitudes and colors of R stars 

The absolute magnitudes My and Mk are calculated as 
follows: 



Mx = mx + 5 - 5 log d - A> 
Ak = 0.114 ^v, 



(where X 



(1) 

VorK)(2) 
(3) 



where d is the distance in pc, and Ay and A-^ are the in- 
terstellar absorptions in the V and K bands, respectively. 
We have adopted the extinction model of Hakkila et al. 
(1997|)]with the ^kMv ratio from |Cardelli et al. {198^ . 
We do not attempt to include circumstellar extinction, 
since this is shown by infrared color to be small or non- 
existent for the R stars. The extinction corrections to the 
K magnitudes are small, less than 0.2 mag. 

Fig. H shows the color - absolute magnitude diagram, 
Mkq versus {V~K)q, for the stars with available K magni- 
tudes. The error bars on Mk are calculated from the upper 
and lower bounds on the parallax listed in Tables & As 
noted in Paper I, the parallax uncertainties are computed 
so as to produce symmetric error bars on the magnitude. 
The resulting uncertainties in the absolute magnitudes are 
often larger than 1 mag (see Fig. ||), and dwarf the uncer- 
tainties on the measured apparent magnitudes and on the 
interstellar extinction. 



Fig. 3. Dereddened absolute magnitude Mkq versus {V — 
K)f) color for R stars. Open symbols: stars of spectral type 
RO - R3; filled symbols: stars of type R4 - R8; squares: 
stars of unknown R subtype. The symbol size is inversely 
proportional to the error on Mk, as calculated from the 
upper and lower bounds on the parallax listed in Tables |^ 
& ||. Large crosses denote stars with H — K > 0.3 that are 
more likely to be N-type rather than R-type carbon stars 
(see Fig. |^). Only stars with K < 8 have been included. 
The ellipse indicates the average location of clump giants, 
according to lAlves (2000| ,see also Fig. |) 



Fig. 1^ shows that the difference between late- and 
early-type R stars, assessed from their photometric, spec- 
tral, variability and mass-loss properties as discussed in 
Sect. lA^ is also apparent in their luminosities: late-type 
R stars with {V — A')o > 4 tend to be more luminous than 
early-type R stars. Late-type R stars thus clearly mark 
the bottom of the asymptotic giant branch. 



3.1. Average absolute magnitude of early-type R stars 

To derive an average absolute magnitude from a set of 
measured parallaxes is an operation plagued by many bi- 
ases (e.g., Arenou & Luri 1999), and several simulations 



have been performed to identify their importance for the 
present sample. 

More precisely, the simulations described below eval- 
uate whether the reduction scheme used in the present 
work is able to correctly retrieve the absolute magnitude 
of the sample stars. In a first step, it is assumed that all 
the 59 stars from the present sample with a K magni- 
tude available have the same absolute magnitude, namely 
M^ — —2.0. In that framework, a true parallax may be 
assigned to each star of the sample, depending on its ob- 
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served apparent K magnitude according to the relation 
log zj = 0.2 (M^ — K) — 1. In the first simulation, the orig- 
inal abscissa residuals of the IAD file were replaced by 
residuals drawn from a normal distribution of zero mean 
and ay = 1 mas (standard deviation of the residuals). 
To artificially increase the sample size, 32 different sets 
of residuals are drawn for the IAD file corresponding to 
each of the 59 R stars, totaling to 1888 simulated stars. In 
this test, the IAD file only serves to define the time sam- 
pling of the observations and the scanning angle on the 
sky f or ea ch observation. The reduction scheme outlined in 
Sect. 2J is then applied to the total sets of modified IAD 
files. Results are shown in the upper left panel of Fig. ^. 
In a second simulation, the reduction scheme is applied 
to the actual abscissa residuals which do not necessarily 
follow a normal distribution with a standard deviation of 
(Tt, = 1 mas. Any differences between the results of the 
first and second simulation will thus highlight the impact 
of the measurement errors. 

The results are shown in the upper panel of Fig. ^, 
and the first three moments of the retrieved distributions 
are listed in Table |[ where the coefficient of skewness 
is defined as /is/cr"^ (and /ia is the centered moment of 
order 3). The method is able to retrieve the true abso- 
lute magnitude when the uncertainties on the residuals 
have a standard deviation of 1 mas. If the actual residuals 
are used instead, the computed distribution is retrieved 
as well, except for a tail extending towards fainter AfK, 
as apparent on Fig. |^. That tail is not visible in Fig. ^, 
though, since stars with K > 8 have very large and noisy 
abscissa residuals and were not displayed. 

Other simulations were carried out as above, except 
that Gaussian distributions with ctq ~ 0.5 and 1 mag 
were adopted for the true absolute magnitude distribution. 
Basically the same features as for the case ao = emerge, 
namely the fact that the retrieved distribution using the 
actual IAD residuals is skewed towards fainter absolute 
magnitudes, with a bias of 0.55 mag for the ao — 1 mag 
case, which however reduces to 0.16 mag if only stars with 
X < 8 are kept (Table |). 

The distribution derived from the true distribution 
{Mf^,ao) = (—2.0,1.0) provides a fair match to the ob- 
served distribution (bottom right panel of Fig. ||), charac- 
terized by a weighted-average absolute magnitude 



-6 -4 



-2 



2 -t 



-4 -2 



Mk 



-1.54, 



(4) 



where MK,i is the individual absolute magnitude of R stars 
with {V - K)o < 4 and fsT < 8, = l/e{MK,if and 
e(MK,i) = 5 (log rot/ - lognjo) = 5 (logtuo - logWi), and 
cc7o, tul, wu as derived by our reduction scheme are given 
in Tables | & |. 

Thus it may be concluded that the early-type R 
star s ha ve a true absolute-magnitude distribution close 
to (M)^,cro) = (—2.0,1.0). Since the present sample of 
R stars is basically magnitude-limited, it is subject to 
the Malmquist bi as favoring high-luminosity stars. The 
Malmquist (1936) bias is probably responsible for the 
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Fig. 4. The distribution of absolute magnitudes Mk for 

the sample of R stars with K < 8 and {V — K)o < 4 ob- 
tained by reprocessing the IAD as indicated in Sect. 2.5 
(thick line) , for different assumptions of the true Mk dis- 
tribution (thin line): upper panel: Af^ = —2.0 for all stars; 
(ii) middle panel: Gaussian distribution of mean M^ = 
—2.0 and (Tq = 0.5 mag; (iii) bottom panel: Gaussian dis- 
tribution of mean = -2.0 and ctq = 1-0 mag. The 
left panels correspond to the simulations performed on 
the IAD files with abscissa residuals drawn from a nor- 
mal distribution of mean zero and (T„ = 1 mas, while the 
right panels correspond to the simulations performed on 
the IAD files with the actual abscissa residuals. In the bot- 
tom right panel, the actual observed distribution (scaled 
by a factor of 30) is represented by the dashed line. The 
bin size in that panel has been enlarged due to the smaller 
size of the observed star sample 



high-luminosity tail and the deficit of low-luminosity stars 
in the observed sample (bottom-right panel of Fig. ||), as 
compared to the simulated Gaussian population. 



3.2. Early R stars as members of the red clump 

Fig. 1^ reproduces the color-magnitude diagram from Fig. |^ 
for R stars with {V ~ K)o < 4 and K < 8, and compares 



them to 238 nearby red clump stars as compiled by Alvei 
|(2000| ). The error bars have been excluded for clarity. 
Red clump stars have become of renewed interest recently 
because Hipparcos has made possible the calibration of 
t heir absolute magnitudes v i a observations of nea rby stars 
( Paczyhski fc Stanek 1998 ; Girardi et al. 1998 ), making 



them a powerful tool for distance measurements within 
the Galaxy and the Local Group. Physically, these objects 
are the analogues of the red horizontal branch stars seen in 
metal-rich globular clusters: helium core burning post-red- 
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Table 3. Weighted average, stan dard deviation and coefficient of skewness of the retrieved Mk distributions for 
different choices of the 'true' (M^,o-o) distributions and input IAD files (see text) 

a. All early R stars: 
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Fig. 5. Mk vs. V — K for early R stars with K < 8 
(filled symb ols), and red clump giants (squares) from 
Alves (20001 ). The ellipse is the same as in Fig. |3| 



giant branch stars. Fig. ^ very strongly suggests that the 
R stars are members of the red clump, since the inferred 
true absolute magnitude M]^ = —2.0 (with an estimated 
intrinsic dispersion of 1.0 mag) of R stars (Sect. 3.1) is 



similar to the average Mk = —1.61 ± 0.03 derived from 
the sample of Alves (2000 ). The V — K colors of R stars 
show a somewhat broader spread (but this may be due to 
contamination by a few N stars), and they are on average 
about 0.5 mag redder than the local clump stars, which 



have < V — K 2.35; this is likely to be due to the pres- 
ence of weak to moderate C2 and CN band absorption in 
the blue region of the spectrum. 

4. The space density of early R stars 

In this section, the space density of R stars will be es- 
timated and compared to that of red clump stars. To 
evaluate the completeness of the sample of early R stars, 
Fig. H displays the cumulative frequency distribution of 
early R stars as a function of distance. To circumvent the 
bias on the derived absolute magnitude distribution (see 
Fig. 1^), the distances have been estimated by assigning 
the absolute magnitude = —2.0 to all early R stars. 
Under the assumption of a uniform density no in the plane 
and an exponential distribution perpendicular to the plane 
with scale height zq, the cumulative frequency distribution 
varies with distance d as 



log N{d) = 4 7mo d^x 



Izq 
2 d 



d_ 

zo 



(5) 



The scale height and the space density of R stars can 
then be estimated from the values that must be assigned 
to the parameters zq and no in order to fit the observed 
distribution. It is found to be about uq ~ 4.5 10~* pc~'^ 
for zq = 300 pc. Fig. ^ reveals that, for this parameter 
set, early R stars follow the predicted uniform/exponential 
distribution up to about 600 pc (involving about 60% of 
the total sample). It needs not be emphasized that these 
values are very uncertain though, as zq is only poorly con- 
strained by the fit and both zq and hq depend on the quite 
uncertain assignment. Adopting M'^ = — 1 instead 
yields no ~ 1.5 10~^ pc~^ and zq ^ 180 pc. 
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Fig. 6. Cumulative frequency distribution of the distances 
of early R stars (thick line), as compared to predictions 
for a stellar sample uniformly distributed in the plane 
and with an exponential distribution perpendicular to 
the plane (thin line). A scale height of 300 pc and a 
true absolute magnitude — —2.0 have been adopted. 
The thin dashed line corresponds to the predictions for a 
spherically-Sjfmmetric homogeneous sample 



We have shown that the early- and late-type R stars 
have different properties. Late-type R stars generally have 
subtypes later than R4, although a better definition uses 
their color indices {V - K)o > A or (H - K) > 0.3. Most 
of these stars are long-period variables of the SR or Mira 
types, and they often exhibit excess emission at 12 fiia due 
to dust, indicative of mass loss. In the HR diagram, these 
stars mark the bottom of the AGB, and are thus related, 
or even identical, to the N stars. 

The early-type R stars form a rather homogeneous 
group, with most of its members populating the red 
clump. They are well represented by a population hav- 
ing a Gaussian distribution of absolute magnitudes such 
that (M^,o'o) = (— 2.,1.). They are not variable and do 
not appear to suffer from substantial mass loss, as judged 
from the absence of 12 /im excess. They represent fewer 
than 0.1% of the clump stars. Their nature as carbon stars 
is not due to low metallicity: R stars have m etal abun- 
dances which are close to solar ( Dominy 1984 ) , except for 
the carbon excess. 

The location of the R stars in the red clump provides 
strong suppor t for the formation scenario suggested by 
Dominy (1984 ): the helium core flash mixes carbon to the 
surface of the star. This mixing episode may perhaps be 
rotationally-induced as the star is spun up by the accretion 
of a former companion. This scenario has the advantage 
of accou nting as well for the lack of binary systems among 
R stars ( |McClure 1997b| ), but it requires a braking mech- 
anism (e.g. magnetic braking as advocated by Habets & 



Zwaan (1989| )?) since R stars do not currently appear to 



The sample of red clump stars by Alves (2000| ) con- 
tains 284 stars with 5% parallaxes from the Hipparcos 
catalogue and is roughly complete to 12 mas or about 
83 pc. The volume density of the red clump stars is thus 
about 1.14 10"'' pc"-^ and the R stars form about 0.04% 
(or 0.13% if = -1) of the total. There are additional 
uncertainties because, for a given star, there is no way 
to tell whether it is a clump star or is on the red giant 
branch. Nevertheless, the qualitative conclusion is clear: a 
very small fraction of clump stars in the solar neighbor- 
hood is carbon-rich. 



5. Conclusions 

In this paper, we have re- analyzed the Hipparcos Inter- 
mediate Astrometric Data for a sample of 83 carbon stars 
of class R, the class of warm evolved carbon red giants. 
A special reprocessing scheme has been used, fitting the 
distance modulus instead of the parallax (thus ensuring 
its positiveness) and rejecting outliers within a bootstrap 
scheme ensuring consistency between the solutions ob- 
tained from the NDAC and FAST data. This reduction 
scheme is not free from biases: Monte- Carlo simulations 
have shown that it tends to produce parallaxes that are 
too large when e{w)/vj ^ 1. Despite these biases, the sim- 
ulations allow the identification of the true properties of 
the population. 



be rapid rotators. The fraction of stars which become R 
stars is much smaller than the fraction (10% - 50%, see 



Wallerstein & Knapp 1998) of stars which become carbon 



stars on the AGB, and it is unlikely to have any significant 
influence on the number of carbon AGB stars. It suggests 
that a small fraction of the RR Lyrae stars and blue hor- 
izontal branch stars in the Galactic halo are also carbon 
rich ~ these stars however are too hot to have C2 in the 
atmosphere. 
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